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Introduction. 

Theoretical  models  for  chemical  lasers  depend  on  a 
variety  of  assumptions  and  empirical  data  to  provide 
closure  and  simplify  solution  of  the  governing 
equations.  Among  the  various  assumptions  and 
empirical  data  that  have  been  built  into  models  for 
chemical  lasers  are  assumptions  regarding  flow 
steadiness  in  the  time  domain  and  geometric  similarity 
of  the  spatial  domain.  The  work  discussed  here  is 
directed  toward  elucidating  and  increasing  the 
understanding  of  these  assumptions  commonly  used  in 
chemical  laser  simulation  and  the  impact  of  their  usage 
upon  the  predictions  of  these  models.  These  efforts  in 
turn  are  directly  linked  to  efforts  to  achieve  improved 
chemical  laser  efficiency  and  performance,  as 
excursions  outside  the  assumed  to  be  ‘well  understood’ 
traditional  operational  parameter  space  are  increasingly 
necessary. 

Problem  and  Methodology. 

Previous  work  by  Madden  and  Miller^  presented 
computational  data  of  the  chemical  oxygen-iodine  laser 
(COIL)  indicating  that  the  flowfield  should  be  marked 
by  substantial  flow  unsteadiness.  Although  COIL 
experiment  data  does  not  exist  to  confirm  this 
prediction,  non-reacting  flow  experiments  at  similar 
flow  conditions  do  substantiate  this  prediction. 
Building  on  the  earlier  results,  the  objective  of  this  work 
is  to  further  elucidate  the  flow  unsteadiness,  its’ 
underlying  physics,  and  the  manifested  effects  in  COIL 
performance. 

The  objective  of  this  work  is  to  generate 
understanding  of  the  temporal  and  spatial  dependence 
of  the  physical  processes  underlying  chemical  lasers 
through  systematic  testing  of  the  models  employed  in 
simulation  of  these  devices.  This  objective  is  achieved 
through  3-D  time  accurate  numerical  simulation  of  the 
COIL  flowfield  and  related  flowfields  from  similar  fluid 
dynamic  experiments. 


The  gas  flows  within  chemical  lasers  can  best  be 
described  as  the  flow  of  particles  of  different  chemical 
composition  with  collisional  interactions  occurring 
between  the  particles  and  between  the  particles  and  the 
photons  within  the  radiation  field.  Mathematically,  this 
flow  of  particles  is  treated  as  a  continuum  and  is 
approximated  by  Navier- Stokes  continuity  equations  for 
mass,  momentum,  and  energy.  In  integral  form,  these 
equations  are  given  by: 

f  III  (2)  +  hdA  =  III  (5)  dvol  [  1  ] 

vol  A  vol 

and  the  hyper  vector  of  the  conserved  variables  (O)  is; 
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The  inviscid  flux  tensor,  F ,  the  viscous  flux  tensor, 
F^,  and  the  hyper  vector  (“S')  of  sources  for  the 
production  and  loss  of  the  reacting  species  are: 
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The  production  rates  due  to  chemical  reaction  of  the 
species  pi  are  determined  for  each  species  continuity 
equation  source  by: 
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7=1 


N 


[6] 


where  the  forward  production  rates  of  each  reaction  are 
determined  from  the  Arrhenius  rate  law  and  the 
backward  rates  determined  from  the  equilibrium  and 
forward  rates: 
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Accurate  description  of  molecular  diffusion  of 
chemically  reacting  species  is  important  in  the  low 
pressure  flowfields  of  chemical  lasers.  The  Ramshaw- 
Dukowicz^  approximation  for  multi-component 
diffusive  transport  is  used  as  opposed  to  higher  order 
approximations  due  to  the  minimal  loss  in  accuracy  and 
lower  computational  cost  with  respect  to  the  higher 
order  methods.  These  equations,  including  both 
concentration  and  pressure  gradient  contributions  are 
given  by: 


Ppi  =  -nniAm 
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Finally,  the  key  laser  quantity  of  gain  is  determined 
from  the  flowfield  solution  using  the  equations: 
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The  computational  fluid  dynamics  (CFD)  code 
GASP  is  employed  to  solve  these  equations  for  the 
simulations  performed  here.  GASP  employs  domain 
decomposition  for  distribution  of  the  computation 
amongst  multiple  processors  on  parallel  architectures, 
with  inner-iteration  methods  used  to  maintain  fiilly 
implicit,  time  accurate  integration  of  the  solutions. 

Reacting  flow,  i.e.  COIL,  and  non-reacting  flow 
simulations  are  performed  in  this  work.  The  GASP 
COIL  model  utilizes  10  species  mass  conservation 
equations  for  the  chemically  reacting  components  of  the 
COIL  flow  in  addition  to  the  base  conservation 
equations  for  momentum  and  energy.  An  effective 
binary  diffusion  model  is  used  to  describe  concentration 
and  pressure  contributions  to  mass  diffusion,  an 
important  process  in  the  low  density  COIL  flowfield.  A 
10  species,  22  reaction  finite-rate  chemistry 
mechanism^  is  used  to  model  the  gas  phase  chemical 
kinetic  processes  that  generate  the  population  inversion 
in  atomic  iodine  necessary  for  laser  oscillation  in  COIL. 
This  model  tracks  the  continuity  of  He,  02(^A),  02(^A), 
02(‘E),  H2O,  I2,  I2*,  l"Pl/2,  l'P3/2,  and  CI2,  a  subset  of  the 
overall  set  of  species  in  the  flow  that  in  conjunction 
with  the  finite-rate  chemistry  mechanism  has  been 
shown  to  accurately  reproduce  the  COIL  chemical 
processes. 

The  computational  grid  used  for  this  first  COIL 
simulation  consists  of  3  blocks  and  3  million  grid  cells. 
The  computational  domain  which  this  grid  discretizes 
represents  the  smallest  geometrically  similar  element 
within  the  COIL  experiment  flowfield  hardware, 
denoted  a  ‘unit-cell.’  The  unit-cell  consists  of  a 
supersonic  (M-2.2)  converging-diverging  nozzle 
section  with  one-half  of  a  large  and  one  small  injector 
orifice  that  issue  reactants  into  the  primary  flow  through 
this  nozzle.  Figure  1  illustrates  the  unit-cell 
computational  grid  within  the  context  of  a  surface 
rendering  of  the  experiment’s  mixing  nozzle.  The 
orifices  inject  a  sonic  mixture  of  He  and  I2  into  a 
subsonic  primary  flow  composed  of  He,  02(^A),  02(^2), 
H2O,  and  CI2  with  the  complex  3-D  flow  structure 
associated  with  the  jet  issuing  from  the  orifice  providing 
the  mechanism  that  mixes  the  two  flows.  Boundary 
conditions  accomplish  the  unit-cell  approximation 
through  the  enforcement  of  planar  symmetry  at  the 
nozzle  centerline  in  the  vertical  direction  and  at  the  side 
boundaries  in  the  lateral  direction.  No-slip  constant 
temperature  boundary  conditions  are  used  at  the  wetted 
surfaces  of  the  nozzle  and  orifices,  with  the  temperature 
fixed  at  400  K  at  the  orifice  walls  and  300  K  at  the 
nozzle  walls.  The  nozzle  subsonic  inflow  boundary 
condition  fixes  the  total  pressure  and  total  temperature 
and  the  species  fractions  at  constant  values,  while  the 
derivative  of  the  pressure  is  set  to  0.  The  nozzle 
outflow  boundary  condition  sets  the  second  derivative 


2 


AIAA-2005-5390 


of  dependent  variables  to  0  as  is  appropriate  for 
supersonic  flows. 

The  non-reacting  flow  simulations  follow  naturally 
from  the  COIL  simulation  and  the  experiments  of  Fric 
and  Roshko"^  discussed  above.  These  models  simulate 
the  same  COIL  hardware  but  do  not  split  the  minimum 
reproducible  set  of  injector  orifices  with  a  symmetry 
plane.  Instead,  the  minimum  reproducible  set  of  one 
large  and  two  small  injector  orifices  are  modeled  in 
their  entirety  with  the  lateral  boundaries  treated  as 
periodic  boundaries  where  waves  exiting  from  one 
boundary  enter  through  the  opposite  boundary  and  vice 
versa.  The  hypothesis  being  tested  here  is  that  the 
symmetry  plane  at  the  centerline  of  the  large  injector 
orifice  interferes  with  the  interaction  of  the  jet  wake 
vortices  with  the  fluid  within  the  jet,  first  developed  and 
explored  in  earlier  work  within  this  line  of 
investigation.^’^  The  computational  grid  used  in  these 
simulations  is  also  somewhat  different  from  that  used  in 
the  first  COIL  simulation.  First,  the  streamwise  extent 
of  the  computational  domain  ends  with  the  end  of  the 
isentropic  expansion  region  of  the  supersonic  nozzle, 
focusing  computational  effort  on  the  injection  region 
and  transonic  region  of  the  throat.  The  second 
difference  is  the  grid  structure,  with  the  new  grid 
utilizing  a  multi-block,  highly  orthogonal  topology  that 
substantially  increases  the  grid  quality.  Finally,  this 
grid  contains  64  million  grid  cells  as  compared  to  the  3 
million  grid  cells  used  in  the  domain  for  the  previous 
simulation.  Fig.  2  illustrates  this  grid. 

The  boundary  conditions  for  this  simulation  remain 
the  same  as  those  used  in  the  first  COIL  simulation, 
with  the  exception  of  the  lateral  boundaries  as  described 
above  to  facilitate  the  capture  of  lateral  fluctuations  in 
the  flow. 

A  sequence  of  physical  models  is  tested  on  this  grid 
to  determine  the  effect  of  varying  composition  and 
chemical  reaction  upon  the  flowfield.  The  motivation 
for  performing  this  test  is  to  determine  if  physical 
models  utilizing  fewer  species  can  be  used  to 
approximate  the  COIL  fluid  dynamics,  and  avoid  the 
computational  expense  of  carrying  the  10  COIL  species; 
computational  costs  scale  approximately  as  the  number 
of  species  squared.  The  models  tested  include  a  single 
gas  He  flow,  a  He,  O2,  I2  flow  that  mimics  the 
composition  of  the  COIL  flow,  and  the  10  species,  22 
reaction  COIL  model  used  in  the  first  simulation. 

References  providing  additional  details  beyond  the 
scope  of  this  paper  may  be  found  in  Madden  et  al.^’^ 

Results. 

We  begin  by  reviewing  and  updating  previous 
results  from  unsteady  simulation  of  COIL.  The  3-D 
GASP  model  for  the  COIL  flowfield  was  executed  in 
time  accurate  mode  utilizing  order  accurate  time 


integration  with  3^^  order  spatial  accuracy  at  a  physical 
timestep  of  1.0x10'^  sec.  The  computation  was 
advanced  to  a  physical  time  of  0.001  sec,  corresponding 
to  100,000  time  steps  in  the  computation.  The  physical 
time  of  0.001  sec  provides  sufficient  advancement  of 
the  computation  to  address  the  characteristics  of  the 
predicted  flow  unsteadiness. 

The  time  accurate  execution  of  the  GASP  COIL 
model  generated  a  prediction  of  the  presence  of  flow 
unsteadiness  that  did  not  decay  over  the  0.001  sec 
interval  that  the  computation  was  advanced.  The 
unsteadiness  was  found  to  extend  from  the  jet/primary 
interaction  region  at  the  point  of  the  He/l2  transverse 
injection  and  continuing  downstream  undiminished. 
Fig.  3  illustrates  the  resulting  impact  of  the  flow 
unsteadiness  upon  the  structure  of  the  flow  within  the 
He/l2  jet.  Here  an  isosurface  of  0.001  constant  I2  mole 
fraction,  a  value  approximately  1%  of  the  concentration 
in  the  pure  He/l2,  is  plotted  within  3-D  space.  The 
surface  demonstrates  the  presence  of  regular,  periodic 
structures  associated  with  unsteady  vortex  generation, 
with  the  structures  present  in  the  flow  from  both  the 
large  and  small  orifices.  Since  molecular  diffusion,  i.e. 
mixing,  and  chemical  reaction  are  strongly  correlated 
with  spatial  gradients  of  reactant  concentration,  the 
presence  of  these  unsteady  flow  structures  is  expected 
to  have  considerable  impact  upon  this  model’s 
predictions  of  device  performance.  The  evidence  of  this 
impact  is  shown  in  Fig.  4  in  the  laser  gain,  or  optical 
wave  amplification  potential  per  unit  distance,  a 
quantity  resulting  directly  from  the  COIL  chemistry. 
Each  position  noted  in  this  plot  corresponds  to  an  X,  Y 
position  with  Z  dependence  integrally  averaged  to 
mimic  the  passage  of  photons  through  the  media  in  the 
manner  exhibited  by  the  laser  itself  The  initial  gain  at 
time  =  0  is  from  steady  state  conditions,  and  as  is 
evident  the  development  of  the  flow  unsteadiness 
substantially  changes  the  magnitude  of  the  laser  gain. 
At  Positions  1  and  2  the  gain  decreases  substantially, 
while  at  Position  3  it  increases  substantially.  The 
explanation  for  this  change  is  the  change  in  the  mixing 
characteristics  and  its’  subsequent  impact  upon  the 
progression  of  the  chemistry.  Fig.  5  shows  the  time,  Z 
direction  integrated  average  of  the  laser  gain  at  various 
points  along  the  X  axis  at  the  nozzle  centerline  in 
comparison  with  measured  laser  gain  and  that  from 
previous  steady-state  simulations.  In  this  comparison 
what  arises  is  a  spread  of  data  due  to  the  influence  of 
advantageous  and  deleterious  processes.  Simulations 
using  reduced  chemistry  mechanisms,  which 
underestimate  the  influence  of  performance  decreasing 
chemical  reactions,  tend  to  align  with  the  measured 
gain.  When  the  deleterious  chemical  reactions  are 
included  as  they  are  in  the  various  full  chemistry 
mechanisms  shown,  the  laser  gain  predicted  falls  below 
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the  measured  gain.  When  the  effects  of  flow 
unsteadiness  are  included  using  the  reduced  chemistry 
mechanism,  the  laser  gain  is  over-predicted  with  respect 
to  the  measured  laser  gain.  A  logical  interpretation  of 
this  plot  is  that  the  full  complexity  of  the  COIL  physics 
is  comprised  of  both  the  full  set  of  chemical  reactions 
and  the  flow  unsteadiness. 

The  impact  of  the  flow  unsteadiness  on  the  optical 
characteristics  of  the  COIL  media  is  a  natural  concern, 
given  the  perturbations  to  the  media  caused  by  the  flow 
unsteadiness.  Fig.  6  shows  the  index  of  refraction,  n, 
variation  within  a  2-D  cut  from  the  COIL  simulation 
perpendicular  to  the  direction  of  optical  wave 
propagation  in  the  device.  A  value  of  I  represents  a 
vacuum,  where  values  greater  than  I  occur  as  optical 
waves  pass  through  matter.  The  greater  the  variation 
and  the  greater  the  value  of  n,  the  larger  the  phase  error 
incurred  in  the  optical  wave  as  it  passes  through  the  gas. 
As  is  shown,  there  are  local  increases  near  the  shock 
waves  in  the  supersonic  region  of  the  nozzle.  However, 
the  values  of  n  remain  very  close  to  I,  with  typical 
values  in  the  5  to  10  cm  region,  where  the  laser  optical 
resonator  resides,  differing  from  I  by  only  I  to  2  xIO'^, 
a  very  small  difference.  This  is  due  to  both  the  very 
low  density  of  the  gas,  and  the  dominant  composition  of 
He  which  has  a  very  low  index  of  refraction  relative  to  a 
vacuum. 

The  work  of  Fric  and  Roshko  indicated  that  a 
variety  of  structures  in  the  jet-in-crossflow  were 
associated  with  the  presence  of  flow  unsteadiness. 
Among  these  are  ‘wake  vortices’  that  project  from  the 
wall  boundary  layer  behind  the  jet.  Three-dimensional 
perspective  can  help  to  further  elucidate  the  wake 
structure  prediction  by  the  simulation.  To  elucidate  the 
flow  structure,  3-D  stream  traces  are  initiated  upstream 
of  the  large  injector  orifice  in  the  XZ  plane  adjacent  to 
the  wall,  as  shown  in  Fig.  7.  In  addition  to  the  stream 
traces,  vortex  cores  in  this  region  of  the  computational 
domain  are  extracted  and  plotted  in  Fig.  7  using  the 
vortex  core  identification  functions  of  Sujudi  and 
Haimes^  as  implemented  in  the  CFD  Analyzer  software 
from  Amtec  Engineering.  While  vortex  cores  are 
captured  and  streamlines  are  shown  projecting 
downward  from  the  wall  toward  the  backside  of  the  jet, 
no  vortex  cores  are  associated  with  this  downward 
projecting  fluid,  indicating  lack  of  rotation. 
Additionally,  this  downward  projecting  fluid  is  adjacent 
to  the  symmetry  plane  at  the  centerline  of  the  large 
injector  orifice,  suggesting  that  the  symmetry  plane  may 
interfere  with  the  development  of  the  ‘wake’  vortex 
structures.  The  absence  of  these  structures  suggests  that 
the  flow  physics  are  not  completely  captured  in  this 
simulation,  and  that  the  captured  characteristics  of  the 
flow  unsteadiness  in  the  simulation  may  not  be  correct. 
Fig.  8  supports  this  notion,  plotting  the  velocity 


component  normal  to  the  symmetry  plane  at  a  variety  of 
monitoring  points  around  the  jet  at  the  first  cell  center 
off  of  this  boundary.  As  can  be  seen,  these  velocities  at 
between  0  and  3  m/sec  are  negligible  in  comparison  to 
the  primary  streamwise  velocity  component  at  order 
200  m/sec  and  the  jet  velocity  of  order  1000  m/sec. 

The  second  simulation  tests  the  hypothesis 
regarding  the  interference  of  the  symmetry  plane  at  the 
midplane  of  the  large  injector  orifice  by  removing  this 
boundary  and  placing  periodic  boundaries  at  the  plane 
midway  between  the  small  injector  orifices,  thus 
encompassing  the  whole  large  injector  orifice  and  both 
small  injector  orifices.  Also,  the  flow  is  composed  of 
the  COIL  diluent  gas  He  alone  to  reduce  the 
computational  cost  for  this  test  that  is  for  examination 
of  the  fluid  dynamic  aspects  of  the  COIL  flow. 

The  GASP  model  for  this  simulation  was  executed 
using  order  temporal  and  3^^  order  spatial  accuracy. 
The  physical  timestep  used  was  I.OxIO'^  sec  for  55,000 
timesteps,  advancing  the  flow  in  time  to  0.00055  sec. 
The  flow  features  captured  by  this  simulation  are  shown 
in  different  perspectives  in  Figs.  9  and  10.  Figs.  9  and 
10  show  a  composite  of  streamtraces  initiated  in  the 
boundary  layer  upstream  of  the  large  injector  orifice, 
vortex  cores,  and  a  Mach  number  =  I  isosurface  in  the 
nearfield  of  the  orifices.  The  streamlines  trace  the 
primary  boundary  layer  flow  as  it  passes  around  the 
large  injector  jet  and  wraps  around  it.  As  the 
streamlines  wrap  around  the  back  side  of  the  large 
injector  jet,  they  project  downward  with  rotation  in  a 
manner  consistent  the  wake  vortex  structures  observed 
by  Fric  and  Roshko.  Vortex  cores  are  extracted  at  the 
center  of  these  structures,  confirming  their 
interpretation  as  wake  vortices.  Comparing  with  this 
same  plot  for  the  reacting  COIL  simulation  in  Fig.  7 
where  vortex  cores  were  not  found  in  conjunction  with 
the  downward  projecting  streamlines,  the  conclusion  to 
be  drawn  is  that  the  symmetry  boundary  condition  at  the 
centerline  plane  of  the  large  injector  orifice  interfered 
with  the  formation  of  the  wake  vortices.  Fig.  10  shows 
the  same  perspective  in  Fig.  9  from  a  greater  distance  to 
provide  greater  detail  regarding  the  interaction  between 
the  wake  vortices  and  the  fluid  from  the  injector 
orifices.  As  the  wake  vortices  project  downward,  they 
impact  the  back  side  of  the  fluid  from  the  injector 
orifice  and  change  trajectory  to  follow  that  of  the  fluid 
from  the  orifices.  This  behavior  is  noted  not  just  with 
the  interaction  of  the  wake  vortices  behind  the  large 
injector  orifice  with  its’  fluid,  but  also  with  those  behind 
the  small  injector  orifices  and  the  small  injector  jet. 
Examining  the  Mach  number  =I  isosurface  in  Figs.  9 
and  10  shows  that  one  set  of  characteristic  vortices 
identified  by  Fric  and  Roshko  in  their  experiments,  the 
‘ring’  vortices  that  initiate  on  the  upstream  side  of  the 
jet  in  the  region  where  the  primary  and  secondary  flow 
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boundary  layer  interact  and  wrap  circumferentially 
about  the  jet,  are  absent  from  this  flow.  The  absence  of 
these  structures  may  indicate  a  fundamental  difference 
between  the  strongly  compressible  COIL  flowfield  and 
the  low  Mach  number,  incompressible  flowfield  in  the 
Fric  and  Roshko  experiment,  with  the  shocks  present  in 
the  COIL  jet  structure  inhibiting  the  formation  of  these 
‘ring’  vortices. 

An  analysis  of  the  fluctuations  associated  with  the 
interaction  between  the  wake  vortices  and  the  jet  is 
shown  in  Figs.  11  -  16  plotting  the  temporal  variation  of 
the  flow  velocity  components  and  fourier  analysis  of 
these  variations.  These  plots  show  that  of  the  10 
monitoring  points  located  in  the  proximity  of  the  jet. 
Points  3  and  4  demonstrate  the  largest  amplitude 
fluctuations.  Both  of  these  points  are  located  on  the 
backside  of  the  large  injector  orifice  jet  where  the 
interactions  with  the  wake  vortices  occur.  As  can  be 
seen,  the  fluctuations  occur  with  all  3  velocity 
components.  Of  special  note  is  the  fluctuation  in  the  Z 
direction  that  is  perpendicular  to  the  symmetry  plane  of 
the  large  injector  orifice.  The  large  amplitude  of  this 
fluctuation,  a  factor  of  10  greater  that  those  in  Fig.  8 
from  the  COIL  simulation  where  the  symmetry 
boundary  is  present,  confirms  the  expectation  that 
removal  of  the  symmetry  boundary  condition  at  the 
centerline  plane  of  the  large  injector  orifice  would 
facilitate  its’  capture  in  the  simulation. 

Fourier  analysis  of  these  fluctuations  indicates 
dominant  frequencies  of  74.8  kHz  and  147.4  kHz  in  the 
U  and  V  velocity  components  for  Points  3  and  4  on  the 
downstream  side  of  the  large  injector  orifice  jet;  Points 
6-9  near  the  wall  behind  the  small  injector  orifice  jets 
have  dominant  frequencies  at  37.4  and  74.8  kHz  in  the 
U  and  V  velocity  components.  The  remaining  points 
follow  the  frequency  distribution  of  Points  3  and  4.  The 
W  or  lateral  velocity  component  of  particular  interest 
here  has  dominant  frequents  of  37.4,  110.0,  and  184.8 
kHz.  Points  6-9  have  these  same  dominant  frequencies 
plus  dominant  frequencies  at  74.8  and  147.4  kHz.  The 
behavior  suggests  that  the  frequencies  associated  with 
the  large  injector  orifice  jet  dominate  the  flowfield,  but 
the  small  jet  orifices  have  their  own  frequencies 
superimposed  upon  the  wave  pattern  generated  by  the 
large  jet.  Since  diameter  of  the  small  injector  orifice  is 
exactly  Z2  of  the  diameter  of  the  large  injector  orifice, 
one  would  expect  resonances  at  the  frequencies  of  the 
large  jet  and  self  resonances  at  frequencies  differing  by 
a  factor  of  2. 

A  separate  execution  of  this  same  model  using  the 
same  execution  configuration  was  performed  using 
factor  of  2  decreases  in  the  number  of  grid  cells  in  each 
spatial  direction,  giving  a  total  of  8  million  grid  cells  as 
opposed  to  64  million  grid  cells  in  the  previous 
simulation.  This  allows  estimates  for  sensitivity  of  the 


calculation  to  grid  density  to  be  performed.  Figs.  17,  18 
show  the  temporal  variation  of  the  W  velocity 
component  and  the  Fourier  analysis  of  this  data. 
Comparing  the  amplitude  of  the  fluctuations  at  Points  3 
and  4  to  the  dense  grid,  we  see  amplitudes  of  18.5  and 
22.5  m/sec  respectively  in  this  grid  versus  21.2  and  27.0 
m/sec  in  the  dense  grid.  The  dominant  frequencies  at 
Points  3  and  4  are  36.0,  108.4,  and  180.4  kHz  versus 
37.4,  110.0,  and  184.8  kHz  on  the  64  million  cell  grid. 
The  better  frequency  resolution  on  the  smaller  grid 
afforded  by  the  longer  physical  time  advancement, 
0.0021  sec  versus  0.0005  sec  on  the  64  million  cell  grid 
did  not  substantially  change  the  frequency  response  of 
the  flow,  nor  did  the  differences  in  spatial  resolution. 
With  regard  to  fluctuation  amplitude  differences 
between  the  grids,  increased  numerical  dissipation  on 
the  8  million  cell  grid  may  account  for  some  or  all  of 
the  differences.  However,  it  is  also  possible  that  small 
shifts  in  the  locations  of  various  fluid  dynamic 
structures  underlying  the  unsteadiness  occur  when 
going  from  one  grid  to  the  other,  changing  the  response 
characteristics  at  the  monitoring  points  that  are  located 
at  the  same  positions  for  both  simulations.  It  should  be 
noted  that  this  same  behavior  of  increasing  fluctuation 
amplitude  with  decreasing  grid  resolution  while 
maintaining  similar  fluctuation  frequency  was  observed 
previously  by  Madden  and  Miller  for  2-D  simulations 
of  the  cylinder-in-crossflow. 

To  test  the  effect  of  flow  composition,  a  simulation 
using  He,  O2,  and  I2  to  mimic  the  COIL  composition 
was  executed  using  the  same  execution  parameters  from 
the  previous  simulations  with  advancement  to  0.00013 
sec.  Figs.  19  -  24  show  the  velocity  component 
temporal  variation  and  fourier  analysis  using  the  64 
million  grid  cell  grid.  The  primary  fluctuation 
frequencies  in  the  U  velocity  component  occur  at  61.5, 
192.3,  and  330.8  kHz;  38.5,  61.5,  84.6,  184.6,  338.5 
kHz  for  the  V  velocity  component;  and  169.2  kHz  for 
the  W  velocity  component  with  Points  3  and  4 
dominating  the  fluctuation  energy.  The  lower 
frequencies  in  the  range  of  10  to  50  kHz  are  not  well 
developed  and  will  be  resolved  by  advancing  the 
simulation  further  in  time.  As  can  be  seen  in  the 
temporal  variation  plots,  the  amplitude  of  the 
fluctuations  is  attenuated  somewhat  in  time,  although 
temporal  advancement  is  not  sufficient  to  determine 
whether  the  amplitudes  will  further  diminish  or  level- 
off  with  time. 

Comparing  the  He-02-l2  flow  simulation  with  the 
He  flow  simulation,  the  100  to  300  kHz  frequency  range 
contains  the  highest  energy  as  opposed  to  the  30  to  80 
kHz  range  in  the  He  simulation,  although  high 
amplitude  components  do  exist  at  110  and  184.8  kHz  in 
the  He  case.  Simulation  of  the  He-02-l2  flow  on  the  8 
million  cell  grid  provides  another  point  of  comparison. 


5 


AIAA-2005-5390 


particularly  with  resolution  of  frequencies  into  the  100 
Hz  range  as  shown  for  the  W  velocity  component  in 
Figs.  25,  26.  As  indicated,  the  most  energetic 
frequencies  captured  on  the  8  million  cell  grid  occur  at 
22.0  and  162.4  kHz  with  the  higher  frequency 
dominating  in  energy  by  a  factor  of  10,  similar  to  the  64 
million  cell  simulation.  Thus  He  flow  simulation  tends 
to  favor  the  lower  frequency  range,  whereas  the  He-02- 
I2  flow  simulation  tends  to  strongly  favor  the  high 
frequency,  attributable  to  differences  in  the  freestream 
velocity  upstream  of  the  jets  which  is  374  m/sec  in  the 
He  flow  and  22 1  m/sec  in  the  He-02-l2  flow. 

Similar  to  the  comparison  of  the  results  on  the  8 
million  and  64  million  grid  cell  simulations  for  the  He 
flow,  the  He-02-l2  flow  simulation  shows  greater 
differences  in  fluctuation  amplitude  than  in  frequency. 
Comparing  the  W  velocity  fluctuations  at  Point  4  shows 
the  amplitude  near  the  end  of  the  physical  time 
advancement  to  be  63.5  m/sec  on  the  64  million  cell 
grid  as  opposed  to  98.5  m/sec  on  the  8  million  cell  grid; 
amplitudes  at  Point  3  were  negligibly  small  for  this 
flow.  Again,  the  effect  of  grid  resolution  on  fluctuation 
amplitude  is  found  to  be  pronounced. 

Finally,  initial  results  for  a  full  reacting  flow  COIL 
simulation  utilizing  a  1  million  grid  cell  version  of  the 
64  and  8  million  cell  grids  discussed  above  are  shown  in 
Figs.  27,  28.  Given  the  observations  regarding  the 
effect  of  decreasing  grid  resolution  on  fluctuation 
amplitude,  the  important  result  is  the  observed 
fluctuation  frequency.  The  peak  fluctuation  frequency 
observed  in  this  simulation  within  the  limits  of  the 
advancement  in  physical  time  is  113.6  kHz,  consistent 
with  the  higher  frequency  fluctuations  found  in  the  He- 
O2-I2  flow  simulation. 

The  attention  paid  to  this  detail,  while  seemingly 
trivial  on  the  surface,  is  justified  based  upon  the  impact 
that  these  fluctuations  have  upon  COIL  performance. 
As  demonstrated  earlier,  the  fluctuations  have  a 
significant  impact  upon  the  laser  gain  through  the 
mixing  and  chemistry.  An  additional  impact  occurs 
through  Doppler  shifts  to  the  frequency  of  radiation  that 
the  atoms  in  the  resonator  field  experience  due  to  the 
fact  that  they  are  translating  with  a  bulk  velocity, 
causing  a  decrease  in  laser  gain,  shown  explicitly  in 
equations  [11]-[13].  Z  direction  fluctuations  are 
significant  for  this  particular  hardware  since  the  primary 
direction  of  transit  of  the  radiation  through  the  device  is 
in  the  Z  direction.  In  similar  COIL  hardware, 
fluctuations  in  the  Y  direction  would  be  significant  in 
this  same  sense  because  the  path  of  optical  transit  is  in 
the  Y  direction. 

Summary  and  Conclusions 

In  conclusion,  the  results  provided  here  establish 
that  the  unit-cell  approximation  as  previously 


implemented  in  simulations  of  the  COIL  flowfield  with 
a  symmetry  plane  boundary  splitting  an  injector  orifice 
interferes  with  the  capture  of  flow  structures  associated 
with  unsteadiness.  Non-reacting  and  reacting 
simulations  of  COIL  hardware  indicate  that  lateral 
passage  of  waves  normal  to  this  symmetry  plane  is 
associated  with  fluctuations  in  the  flow  about  this  plane. 
An  examination  of  the  effect  of  variations  in  flow 
composition  from  He  to  He-02-l2  to  the  full  COIL 
species  set  showed  that  the  He  flow  tended  to  fluctuate 
predominantly  in  the  30  to  80  kHz  range  whereas  the 
He-02-l2  and  COIL  flows  fluctuated  in  the  100  to  300 
kHz  range.  In  the  reacting  flow  COIL  simulation  where 
the  symmetry  plane  was  present,  the  flow  unsteadiness 
found  there  was  shown  to  have  a  substantial  effect  on 
the  laser  gain  prediction. 

Future  work  will  continue  the  line  of  investigation 
initiated  with  the  reacting  flow  COIL  simulation  using 
the  full  unit-cell  as  well  as  larger  domains  incorporating 
adjacent  injector  orifices  laterally  and  across  the  nozzle 
centerline. 
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Figure  1.  Computational  grid  used  in  the  first  3-D 
GASP  COIL  simulation,  utilizing  a  symmetry  plane 
boundary  condition  at  the  center  of  the  large  injector 
orifice.  The  grid  is  shown  in  relation  to  a  rendering  of 
the  wetted  surfaces  in  the  COIL  laser  hardware. 


dynamic  structures  manifested  in  the  He/l2  jet  from  the 
first  3-D  GASP  COIL  simulation  utilizing  a  symmetry 
plane  at  the  center  of  the  large  orifice  (see  Fig.  1).  A 
surface  of  constant  I2  mole  fraction  acts  as  a  scalar  tracer 
for  the  jet  structure. 


Figure  2.  Computational  grid  used  all  other  3-D 
GASP  simulations  of  COIL  hardware,  shown  in  relation 
to  a  rendering  of  the  wetted  surfaces  in  the  COIL  laser 
hardware.  Note  that  the  lateral  boundaries  are  now 
between  the  small  orifices.  Also  note  the  use  of 
multiblock,  highly  orthogonal  grid  structure  to  improve 
grid  quality  over  the  grid  used  in  the  first  simulation. 


Figure  4.  Time  traces  for  the  Z-averaged  laser  gain  at 
different  positions  in  the  laser  resonator  section  of  the 
simulated  COIL  hardware.  This  data  was  taken  at  a 
sampling  rate  of  10  iterations  of  the  3-D  GASP  COIL 
simulation  (see  Fig.  1),  for  a  temporal  sampling  rate  of 
1x10'^  sec. 
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Comparison  of  Simulation  Probe  Beam  Averaged  and  Measured 
Gain  at  the  Nozzle  Centerline. 


x(cm) 

Figure  5.  Comparison  of  time,  Z-averaged  laser  gain 
from  the  first  unsteady  COIL  simulation  (see  Fig.  1)  with 
experiment  data  and  previous  steady-state  COIL 
simulations.  The  steady  state  simulations  include 
varying  levels  of  chemistry  modeling  fidelity. 


Figure  6.  Index  of  refraction  within  the  2-D  plane  at  the 
center  of  the  large  injector  orifice  from  the  first  3-D 
GASP  COIL  simulation  (see  Fig.  1). 


Figure  7.  3-D  streamtraces  (in  black)  initiated  upstream 
of  the  large  injector  orifice  in  the  first  XZ  plane  (0.00098 
cm)  from  the  wall  in  conjunction  with  vortex  cores  (in 
red)  extracted  from  the  first  3-D  GASP  COIL  simulation 
(see  Fig.  1). 


Figure  8.  Lateral  or  Z  velocity  component  (W)  variation 
with  time  at  monitoring  points  located  along  the 
symmetry  plane  in  the  proximity  of  the  injector  orifices 
from  the  first3-D  GASP  COIL  simulation  (see  Fig.  1). 
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Figure  9.  Streamlines,  in  black,  vortex  cores,  in  red,  and 
Mach  number  =1  isosurface  in  the  jet  nearfield  from  the 
He  flow  COIL  hardware  simulation  (see  Fig.  2).  This 
simulation  does  not  use  the  symmetry  boundary  at  the 
center  of  the  large  injector  orifice  as  in  the  first 
simulation,  but  uses  periodic  boundaries  outside  of  the 
small  injector  orifices. 


Figure  11.  Time  variance  of  the  streamwise  or  X 
direction  velocity  (U)  component  at  various  monitoring 
points  in  the  jet  nearfield  from  the  He  flow  COIL 
hardware  simulation  shown  in  Figs.  9  and  10. 


Figure  10.  Streamlines,  in  black,  vortex  cores,  in  red, 
and  Mach  number  =1  isosurface  in  the  jet  near- field 
from  the  He  flow  COIL  hardware  simulation  (see  Fig.  2). 
This  simulation  does  not  use  the  symmetry  boundary  at 
the  center  of  the  large  injector  orifice  as  in  the  first 
simulation,  but  uses  periodic  boundaries  outside  of  the 
small  injector  orifices.  This  perspective  is  further 
removed  than  with  Fig.  9. 


frequency  (Hz) 

Figure  12.  Fourier  analysis  of  temporal  variation  of  the 
streamwise  or  X  direction  velocity  (U)  component  in 
Fig.  11. 
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Figure  13.  Time  variance  of  the  vertical  velocity  (Y) 
component  at  various  monitoring  points  in  the  jet 
nearfield  from  the  He  flow  COIL  hardware  simulation 
shown  in  Figs.  9  and  10. 


component  at  various  monitoring  points  in  the  jet 
nearfield  from  the  He  flow  COIL  hardware  simulation 
shown  in  Figs.  9  and  10. 


frequency  (Hz) 

Figure  14.  Fourier  analysis  of  temporal  variation  of  the 
vertical  velocity  (Y)  component  in  Fig.  13. 
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Figure  16.  Fourier  analysis  of  the  lateral  velocity  (Z) 
component  in  Fig.  15. 
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Figure  17.  Time  variance  of  the  lateral  velocity  (Z) 
component  at  various  monitoring  points  in  the  jet 
nearfield  from  the  8  million  grid  cell  He  flow  COIL 
hardware  simulation. 


points  in  the  jet  nearfield  from  the  He-02-l2  flow  COIL 
hardware  simulation. 
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Figure  18.  Fourier  analysis  of  the  lateral  velocity  (Z) 
component  in  Fig.  17. 
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Figure  20.  Fourier  analysis  of  temporal  variation  of  the 
streamwise  or  X  direction  velocity  (U)  component  in 
Fig.  19. 
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Figure  21.  Time  variance  of  the  vertical  or  Y  direction 
velocity  (V)  component  at  various  monitoring  points  in 
the  jet  nearfield  from  the  He-02-l2  flow  COIL  hardware 
simulation. 


Figure  23.  Time  variance  of  the  lateral  or  Z  direction 
velocity  (W)  component  at  various  monitoring  points  in 
the  jet  nearfield  from  the  He-02-l2  flow  COIL  hardware 
simulation. 
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Figure  22.  Fourier  analysis  of  temporal  variation  of  the 
vertical  velocity  (Y)  component  in  Fig.  21. 
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Figure  24.  Fourier  analysis  of  the  lateral  velocity  (Z) 
component  in  Fig.  23. 
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Figure  25.  Time  variance  of  the  lateral  or  Z  direction 
velocity  (W)  component  at  various  monitoring  points  in 
the  jet  nearfield  from  the  He-02-l2  flow  COIL  hardware 
simulation  on  8  million  grid  cells. 


Figure  27.  Time  variance  of  the  lateral  or  Z  direction 
velocity  (W)  component  at  various  monitoring  points  in 
the  jet  nearfield  from  the  COIL  hardware  reacting  flow 
simulation  on  1  million  grid  cells. 
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Figure  26.  Fourier  analysis  of  the  lateral  velocity  (Z) 
component  in  Fig.  25. 
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Figure  28.  Fourier  analysis  of  the  lateral  velocity  (Z) 
component  in  Fig.  27. 
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